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Motivated by the desire to better understand the class of quantum operations on bipartite systems 
that completely preserve positivity of partial transpose (PPT operations) and its relation to the class 
LOCC (local operations and classical communication), we present some results on deterministic bi- 
partite pure state transformations by PPT operations. Restricting our attention to the case where we 
start with a rank K maximally entangled state, we give a necessary condition for transforming it into a 
given pure state, which we show is also sufficient when K is two and the final state has Schmidt rank 
three. We show that it is sufficient for all K and all final states provided a conjecture about a certain 
family of semidefinite programs is true. We also demonstrate that the phenomenon of catalysis can 
occur under PPT operations and that, unlike LOCC catalysis, a maximally entangled state can be a 
catalyst. Finally, we give a necessary and sufficient condition for the possibility of transforming a rank 
K maximally entangled state to an arbitrary pure state by PPT operations assisted by some maximally 
entangled catalyst. 



I. INTRODUCTION 



In the theory of quantum entanglement the familiar 'distant labs' scenario motivates the restriction of allowed op- 
erations to those that can be carried out by local operations and classical communication (LOCC). Since no entangled 
state can be created from an unentangled one by LOCC, we regard entanglement as a resource in this context. In the 
case of the bipartite pure states, the picture of how this resource can be quantified and how it can be transformed (by 
LOCC) is quite well developed: For bipartite pure states the entanglement of distillation, E d , and the entanglement 
cost, E c , have been shown to both be equal to the entropy of entanglement (for a review see Plenio and Virmani [1]). 
Since Ed and E c are equal, LOCC transformation of bipartite pure state entanglement is asymptotically reversible, 
and the entropy of entanglement is a unique measure of the amount of entanglement available in the many copy 
limit. 

For mixed bipartite states, entanglement transformations are not in general asymptotically reversible, and con- 
sequently, LOCC no longer induces a unique measure of the amount of entanglement (in the asymptotic regime). 
For instance, there are bound entangled states (which by definition have zero entanglement of distillation) which 
have been shown to have non-zero entanglement cost [2]. This has motivated the search for non-trivial extensions of 
LOCC, with respect to which entanglement transformations are asymptotically reversible for all bipartite states. Re- 
cently, Plenio and Brandao [3] proved that the set of operations which asymptotically cannot generate entanglement 
is such a set, with the regularised relative entropy of entanglement as the corresponding measure of entanglement. 
An open conjecture [4] is that the set of PPT operations [5] (defined later in this section) also renders bipartite entan- 
glement transformations asymptotically reversible (in [4] it is shown that it does so for the anti-symmetric Werner 
state). 

For exact, finite transformations, the problem of determining when a given pure bipartite state can be determinis- 
tically converted into another by LOCC is completely solved by Nielsen's majorization theorem [6]: The process is 
possible if and only if the Schmidt coefficient vector of the initial state is majorized by that of the final state. 

In the present paper we investigate exact, finite transformations of pure states by PPT operations. This topic was 
first treated by Ishizaka and Plenio [7]), where the emphasis is on conversion of LOCC-inequivalent forms of multi- 
partite entanglement (e.g. GHZ and W states for three parties) by PPT operations. In the present work we are only 
concerned with bipartite states. 



'Electronic address: william.matthewsObris.ac.uk 
^Electronic address: a. j .winter@bris .ac.uk 



2 



Here, we introduce some conventions that will be used throughout the rest of the paper. Logarithms are always 
taken to base two. If the deterministic transformation of a state p into p' can be accomplished by the class of opera- 

OP OP 

tions OP then we write p — > p'; if it cannot be then we write p -» p' . In discussing the transformation of bipartite 
pure states by any class of operations which contains LOCC, we need only consider the Schmidt coefficients of the 
states since states with the same Schmidt coefficients are equivalent up to local unitary transformations (which are 
obviously contained in LOCC). Since we are only concerned here with classes that include LOCC, we will use the 
state 



PA = EEv / ^7N0'il (!) 

i=l 3 = 1 

(where \ij) = \i)a ® \])b for orthonormal bases b} for Alice and Bob's quotients of the bipartite Hilbert 

space) as a representative of all states with Schmidt coefficient vector A without loss of generality. We use A T (A^) to 
denote the vector obtained by putting the components of A in non-decreasing (non-increasing) order. 

Using this notation, Nielsen's theorem is 

LOCC . , 

P\ -> <^ x -< ( 2 ) 
where the majorization relation is defined on vectors in R d whose components sum to one by 

i 3 
i=l i=l 

We shall use $>k to denote a maximally entangled state of rank K (where we assume that K > 2). 

We will use the symbol X T to denote the partial transpose an operator X. We define the (linear) partial transpose 
map by 

\im\ T = \u){kj\. (4) 

Clearly this is basis dependent, but the eigenvalues (and hence the positivity) of the partial transpose of an oper- 
ator does not depend on this basis choice. For convenience we will use the same basis that we are using as our 
representative Schmidt basis. 

The set of PPT operations on a bipartite system is the set of completely positive trace-preserving (CPTP) maps C 
such that the composition r o C o T is also completely positive, where T denotes the partial transpose map p — > p r . 
Equivalently, PPT maps are CPTP maps that completely preserve the PPT property of states in the same sense that 
completely positive maps completely preserve non-negativity of states: Any extension of a PPT map on a system Q 
onto a larger system QR where we apply the original map on Q and the identity map on R is PPT preserving (it maps 
PPT states to PPT states). 

We make frequent use of the Renyi entropies: For t £ [0, oo] the Renyi entropy at t is defined by 



iVog(Eti^) forte (0,1) U (l,cx.), 



St (A) = <^ lo S l A l for * = °> (5) 

* H{\) for t= 1, 

— log Aj for t = oo, 

where |A| denotes the number of non-zero components of A (i.e. the Schmidt rank) and if (A) = — J2i ^» 1°S ^» ls tne 
Shannon entropy of A. 

In the next section we discuss some necessary conditions for general bipartite pure state transformations by PPT. 
After that we provide part of the theory analogous to Nielsen's theory for deterministic transformations of pure 
bipartite states by PPT operations. In particular, we look at the case where the initial state is maximally entangled. In 
section III we provide a necessary condition for transformations of this type, which we conjecture is also sufficient. 
We make use of this result in section IV to show that the phenomenon of catalysis [8] can occur for PPT operations 
and give a necessary and sufficient condition for when this can occur if both the initial state and catalyst state are 
maximally entangled. We conclude with some suggestions for future work, including a conjecture for extending this 
condition to arbitrary initial and catalyst states. 
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II. GENERAL PPT PURE STATE TRANSFORMATIONS 



Determining whether a particular pure state transformation can be carried out by a PPT operation can in fact 
be formulated as a semidefinite programming problem [7], the difficulty is in phrasing the relevant constraints in 
terms of conditions on the Schmidt coefficients of the pure states. Here we show that certain Renyi entropies of 
the Schmidt coefficient vectors correspond to operationally motivated PPT monotones (allowing us to give some 
necessary conditions for transformations). 

For completeness let us define these operational monotones for any class of operations X: Let p be a density 
operator in B(Ha ®Wb), where B(H) denotes the set of hermitian operators on the Hilbert space H. 



Definition 1. The distillable entanglement of p under the class of operations X is defined by 



and the exact distillable entanglement by 



Ej(p):= sup (lim 10 ^ lim \\* Kj - C j (p^)\\ 1 = o} 



(6) 



2&(p) := lim max 



1*^.-^(^)111 = , (7) 



where in both cases (£j)j £N is a sequence of completely positive trace preserving (CPTP) maps such that each map, Cj, takes 
B ((Ha <8> H B f nj ) to B (C®*' <8> C®*') and belongs to the class X. 

Definition 2. The entanglement cost of a state p, we define by 

E?(p) := iuf { lim lim \\Cj($ Kj ) - p® n > ||i = o) , (8) 



and the exact entanglement cost by 



El(p) := limmin/ 1 ^ 



{Cj^K^-p^W^O , (9) 



where in both cases (Cj) je n is a sequence of CPTP maps such that each map, Cj, takes B (C®-^ 3 <g> <C® K] ) to 
B ((Ha and belongs to the class X. 

Proposition 3. The entanglement cost E^ FT (p\) and distillable entanglement E FPT (p\) of px are both equal to 5i(A), the 
entropy of entanglement of the state. 

Proof. It is clear and well-known that by LOCC operations, both the entanglement cost and distillable entanglement 
of p\ is Si (A). By elementary results of Rains' theory of PPT distillation [5], PPT operations cannot asymptotically 
increase the number of EPR pairs. Hence 

£ c PPT (pa)>S!(A), EY T (px)<S 1 (X). 

Since the opposite inequalities are trivial by LOCC C PPT, we conclude that 

£ PPT (pa) = £ PPT (pa) = Si(A), 

thus coinciding with the corresponding LOCC entanglement cost and distillable entanglement: the entropy of en- 
tanglement. □ 

Proposition 4. The exact distillable entanglement -B ppt (pa) of px is given by S (X) (A). 

Proof. We use a result of Rains [5] on the maximum fidelity obtainable from via PPT operations, to a maximally 
entangled state of Schmidt rank K: by this result, exact transformation (fidelity 1) is possible if there exists an 
operator F with 

PX<F<1, -il <F r < il. 
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Thus, an upper bound to K is given by 

nmx{||F r ||^ \p x <F<t}. 

We claim that this is 1/X[: for assume any F as above, and write it F = p xi + A, A > 0. 

||F r || oo >Tr(|llXll|( p r i+ A r )) 

= A|+Tr(|ll)(ll| r A) (10) 
= Aj+Tr(|ll)(ll|A) > \[. 

Hence H^H^ > X{, with equality achieved for F = p\. 

This gives E^ T (p\) < S 00 (X). The opposite inequality comes from LOCC operations asymptotically achieving 
this bound [9]. ' ' □ 

Proposition 5. The entanglement cost E^ T (px) of p\ is Si/ 2 (X). 

Proof. This is merely an application of the more general result of Audenaert et al. [4] to pure states. □ 

All of these quantities are clearly PPT monotones so they provide a necessary condition for the possibility of a PPT 
pure state transformation: 

PPT 

Proposition 6. If p^ — ► p\ then 

Si/ 2 (ji) > Si/ 2 (A), (11) 
Si(p) > Si (A), (12) 
Soo(p) > Soo(A). (13) 

As Schur concave functions, the Renyi entropies of the Schmidt coefficients at all values of t e [0, oo] are monotones 
for LOCC state transformations. Under PPT however, the Renyi entropies for < t < 1/2 are not monotones. 

Example 7. Consider a pure state p\ with Schmidt spectrum 

._ / 1 1 14 4 9 
"~ V20' 20' 20' 20' 20' 20 

It is easily verified that S 1 / 2 {X) = log 5, and indeed, in accordance with the Proposition 5 the transformation 



^®(n+o(n)) ^ 

5 * A 



is possible by PPT operations for all sufficiently large n. 
But for < a < 1/2, 

log5<5 t (A). 

Since there are trivial examples of initial states such that the opposite inequality is true (e.g. $ 6 ), the Renyi entropies atO <t < 
1 /2 are not PPT monotones. 

S (X) is just the logarithm of the Schmidt rank, so, as was noted in [7], PPT operations can increase the Schmidt rank of pure 
states, a thing LOCC transformations cannot even do with nonzero probability! 

It should be noted that a necessary condition for pure bipartite state transformations by separable operations was 
recently given by Gheorghiou and Griffiths [10]. 



III. PPT TRANSFORMATIONS FROM MAXIMALLY ENTANGLED STATES. 



Unless otherwise stated, the final state is p\, where the Schmidt coefficient vector A e M d is assumed without loss 
of generality to have no vanishing components (so d is the Schmidt rank of the final state). Since any state with 
Schmidt rank not greater than K can be produced from <& K by LOCC, the interesting case for PPT transformations 
is when the Schmidt rank is increased. 
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PPT 

Lemma 8. For any (pure or mixed) final state p, <& K — > P if and only if the solution to the semidefinite program 

min{Tr (P) \P>0,-(K- 1) P r < p T < (K + 1) P r }, (14) 
where P is an hermitian operator on the same Hilbert space as p, is less than or equal to one. 

Proof. The argument is almost the same as the one given in [4], but we present it here for convenience. If there is a 
PPT map C such that C ($k) — P then the map £ made by preceding C with the twirl operation [5] T (which can 
be implemented by LOCC) is also PPT, and does the same transformation since T (<&k) = By symmetry it is 
always possible to write the new map in the form 

£'(r)=FTr($ Jf T) + GTr((l-$jf)r) (15) 

where, in order for the map to be CPTP, F and G must be density operators. In order that £ (r) = p, we require 
F = p, so the desired PPT transformation is possible if and only if there is a density operator G such that the resulting 
£ is PPT. £ is PPT preserving if and only if the operator 

£ (r r ) r = p r Tr (^t) + G r Tr ((1 - <S> K f r) (16) 

= p r Tr ((S - A) t) /K + G r (((1 - l/K) S + (1 + l/K) A) r) (17) 
= Tr (St) (p r /K + (1 - l/K) G r ) + Tr {At) ((1 + l/K) G r - p T /K) (18) 

is positive semidefinite for all positive semidefinite r, where S and A are the projectors onto the symmetric and anti- 
symmetric subspaces of the bipartite space. Since SA = 0, this condition holds if and only if p T j 'K+(l - 1 / 'K) G T >0 
and (1 + l/K) G T — p r / K > 0. We also require G > and TrG = 1. The set of operators which satisfy these 
constraints is precisely the set of feasible points of the SDP (14) which have trace 1. If Pq is feasible, then points 
Pq + tt are also feasible for all t > 0, so if an optimal point has trace not greater than one, it ensures the existence 
of a feasible point which satisfies all the constraints on G (so the transformation is possible); if an optimal point has 
trace greater than one, then clearly no such point exists and hence the transformation is not possible. 

□ 



Proposition 9. 



<5>K Pf ? Px T(K;X)<1, (19) 



where we define T (K;X) to be the solution to the semidefinite program 

id d /A~A~ /XX" 



i>j i>3 



i>j i>j 



with 



a . . = f + \ji))((ij\ + 0"*l)/2 when i + j, 
13 \ when i = j, 

<*H =(\ij)-\ji))({ij\-{ji\)/2, 
and {sij\l < j < i < d}, {atj\l < j < i < d} together constitute d 2 real variables. 
Proof. 



Pa = V x i x j a ij - ^2\/^ x i a ir ( 21 ) 

i>j i>j 



Let II be the projection map on the space of hermitian operators given by 



II (t) = ^ (J i] T(T ij + X] a i3 Ta i3- ( 22 ) 

i>j i>j 
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To show that this map preserves positivity of partial-transpose, we note that 

n (^ r ) r = \ X>'My> (ly'Xy I + b*)0*l) + * (E l"X«l) + O ^E Iw'XjjI j • (23) 

Clearly then, II is positive and preserves positivity of partial transpose. Since p v lies in the image of II, if P is a 

feasible point of the semidefinite program (14) then II {P T ) T is also a feasible point. Since II is also trace reducing, it 

will not change the optimal value of (14) if we impose the additional restriction that P = II (P r ) T ■ This restriction is 
obeyed if and only if we can write P in the form 

p = e s i i °ij + E °« ■ (24) 

i>j i>j 

Substituting this into (14), we obtain (20). □ 
Lemma 10. 

T(K;X)< (2 Sl / a W-l)/(K"-l). (25) 

Proof. The point 

Sij = yfKXj/ {K - 1) , aij = (K - 1) (26) 

is primal feasible since 

E + E + a ^ T = E ^"X"! + E (ly'Xy I + b*)0'*l) > ( 27 ) 

and the other inequalities are obviously satisfied. The primal objective function at this point is 2 Sl ' 2 ^ / (K - 1). □ 

Since the semidefinite program (20) is strictly feasible (take the point = 2, = 2, for example), its solution is 
equal to the solution of the dual SDP [13] so we have 

Proposition 11. 

I _ 1 /„ . A A „ / . \ \ \ 

{l*,t)€R} (28) 



i>j i>j 



R := { GM) 



/'., < 1, < 1, E + E "J"- ( 29 ) 

i>J i>3 



Here {nij | 1 < j < i < d}, {tij | 1 < j < i < d} together constitute d 2 real variables. 
Proof. The semidefinite program dual to (20) is 

maximizc _L_ £ (1 - Tr {Za T i3 )) + ^ E (l - Tr (Za£)) (30) 



i>3 i>3 

subject to Z > (31) 

Tr (Zo£) < 1 (32) 

TV (Zc^) < 1. (33) 

If Z is feasible point of this program then so is II (Z r ) r (where II is the map defined in the proof of Proposition 9). 
So the substitutions 

Z = E % a i3 + E tl 1 "iJ ' X% 3 = 1 ~ 1*0 ' ^ = 1 ~ tl 3 ^ 

i>3 i>3 

result in an SDP with the same solution and this is the one given in the proposition. □ 
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The dual objective at any feasible point of the dual semidefinite program (28) is a lower bound on T (K; A), and 
therefore 

Corollary 12. 

T (K; A) > (k2 s v*W _ / (#2 _ ( 35 ) 

Proof. The point pn 3 — 0,tij = is clearly dual feasible. □ 

Theorem 13. The optimal value of the dual objective that can be attained by a dual feasible point which satisfies the additional 
constraint that 

rank ^ ^ij + Yl *ij a Fj ) = ( 36 ) 

\ i>3 i>3 J 



is given by 



K 2 s i/2W _ i k i , v - ■ . v- 



if 2 -l (if 2 - 1) (K + c* - d) 
where c* is the smallest number c e {1 + d - K, d — 1} satisfying 



v.i=l / i=l 



If none of the integers in the range satisfy this relation then c* = d. 
(The proof of this theorem is given in the appendix.) 



— PPT 

Remark 14. Clearly 7\ (if; A) < T (if; A), so a necessary condition for the transformation <& K — > Pa zs Ti (if; A) < 1. 
Corollary 15. 1/ 04 ^ / (if - 1) > f/zen 



vi=l 



T (if; A) = Ti (X; A) = (V 1 /^) - l) / (if - 1) . (39) 

Proof. In this case c* = d and so, T\ (if; A) = (2 Sl /2( A ) - l) / (if - 1) which, by Lemma 25, is an upper bound on 
T (if; A) . Since Ti is also a lower bound on T, the result follows. □ 

Corollary 16. If Si/ 2 (A) = log if and d > K then the transformation is possible only in the trivial case where the goal state 
is also a maximally entangled state of rank if. 

Proof. The value c = d — 1 satisfies (38) provided that A^ > 1 /if or equivalently if Soo (A) < log if. Since S t (A) is a 
non-increasing function of t, this condition is indeed satisfied. Using this value of c in (37) yields the lower bound 

k( 



which is clearly greater than 1 (implying the impossibility of the transformation) except where (A) = log if which 
(together with S 1 / 2 (A) = log if) implies A is the uniform distribution of size if, so the goal state is a maximally 
entangled state of rank if. □ 

Proposition 17. In the case where the goal state has Schmidt rank three, T (2; A) = Ti (2; A). 
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Proof. To simplify notation we shall here assume that Aj = A}. If \f\i + > y/X^ then we can apply Corollary 15 
and we're done, so we assume that >/Ai + \/X~2 < V^3- In this case 

Ti = (l+4(v/A^+ V^AT) +8v/a7a^) /3. (41) 

We shall show that T is the same by constructing a primal optimal solution. Let 

*11 - \m + gV^A^,«22 = ^V 7 ^,^ = (42) 



'12 



1 4 



2 J 

, , \Al-^2,Si 3 = -x/AiA 3 ,S23 = g V%^3 



(43) 

V^A~ (44) 



S + X] a « a « = p i © p 2 (45) 



p i = 1 E V^KiX^I + | (I23>(23| + |32)(32|) > 



3 

Written as a matrix in the basis, 

: / 3A : + 4 % /A 2 A 3 ~ -3\A^ -3 V / A7A 3 " 

P2 = g I — 3\/A2Al 3A2 + 4V^2A3 — \/A2A3 

\ — 3\/A3Al — \A3A2 3A3 

which can be seen to be positive semidefinite by the Sylvester criterion. Therefore Pi P 2 > 0, and since the other 
primal constraints are clearly satisfied, the point (s*,a*) is primal feasible so 

Ti (2; A) < T (2; A) < £ 4 + ]T a?- = (l + 4 (yX^ + + 8V^) /3 = Ti (2; A) . (46) 

Therefore T (2; A) = Ti (2; A). □ 

It follows directly from this proposition that: 

Theorem 18. A pure state of Schmidt rank three can he produced from <j> 2 (an EPR pair) by PPT operations if and only if its 
Schmidt coefficients obey 

2(v // 44 + \Aw)+ 4 v // ^^ 1 - ( 47 ) 

This fact, and some suggestive numerical evidence, leads us to make the following conjecture 

Conjecture 19. In the dual program (28), the optimal value of the dual objective function is always attained by a point satisfying 
the rank constraint (36) and as a consequence T (K; X) = Ti (K; A). 

IV. CATALYSIS 

Nielsen's theorem was used to show that a phenomenon analogous to chemical catalysis can occur in LOCC 

LOCC 

entanglement transformation [8]. That is, there exist pairs of states p\, p M such that p\ «/» p M but where there is a 

LOCC 

third state such that p\ <g> —> p^® p^. 

A necessary and sufficient condition for the existence of a deterministic catalytic transformation for bipartite pure 
states was recently given in papers by Turgut [11] and Klimesh [12]: 
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(1/3, 1/3, 1/3) 




Accessible Without Catalysis Initial State 



Figure 1: Pure states of Schmidt rank three accessible deterministically from a single EPR pair by PPT and by PPT assisted 
by a maximally entangled catalyst. 



Theorem 20. Given two pure bipartite states p\ and p^, where A T ^ an & \ an ^ ^ & on 't froth have components equal to 



zero, there exists a pure state p^ such that p\ <8> p^ 



LOCC 



Py, ® Pi if and only if the following conditions are satisfied 



S t (X)>S t (p) forte (0,00), 
ft(X)>ft(p) forte (-oo,0] 



where 



jh log (Eti Af) 

— oo 



when t > and A, f= for all i € {1, d}, 
when t = and A, f^ /or a// i e {1, d}, 
otherwise. 



(48) 
(49) 



(50) 



Using the results established in the last section, we can show that catalysis can also occur under PPT operations. 



Theorem 21. $c<&®k 



PPT 



$c ® P\ (where p\ ^ <&x)for some sufficiently large value ofC if and only if Si/ 2 (A) < log K. 



Proof. For a given value of C, the transformation is possible if and only if T (KC; A <g> Uc) < 1/ where Uc denotes the 
uniform probability distribution vector with C elements. As C tends to infinity both the upper bound (25) and lower 
bound (35) on T (KC; A <g> U c ) tend to 2 5l /2(A) jk. Therefore, if S 1/2 (A) < log K then for some sufficiently large C, 
T (KC; \®Uc) < 1 so the corresponding transformation is possible. That the condition is necessary follows from 
Proposition 5 (Si/ 2 cannot increase in any PPT pure state transformation and the Renyi entropies are additive) and 
Corollary 16 (deals with the case where S\/ 2 stays the same). □ 



In the case where the goal state has Schmidt rank three, the states satisfying 



\\\X + ^\[\ 2 ) <i 



T \ T 



(51) 



are exactly those that can be reached from <f> 2 by PPT operations when maximally entangled catalysts of arbitrarily 
high rank are available, according to Theorem 21. This is a strict superset of those Schmidt rank three states which 
can be obtained from <& 2 without a catalyst (see Theorem 18). These regions are illustrated on one cell of the simplex 
of Schmidt coefficient vectors in Fig. 1. 



V. CONCLUSIONS 



We have provided a necessary condition for the exact preparation of a pure bipartite state from a maximally 
entangled state by PPT operations in terms of the Schmidt coefficients of the final state. We conjecture that this 
condition is also sufficient and have shown that this is true when the final state has Schmidt rank three. We have 
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demonstrated that the phenomenon of catalysis occurs in the context of PPT operations. A notable difference from 
LOCC catalysis is that maximally entangled states can act as catalysts under PPT operations (this is impossible with 
LOCC). In the case where both the initial state and the catalyst are maximally entangled, we have given a necessary 
and sufficient condition for the production of a pure state: the Renyi entropy at 1/2 of the Scmidt coefficient vector 
must decrease, unless the final state is the same as the initial state. To give a direction for future work - comparing 
Theorem 21 and the observations on PPT monotones given in Section II with Theorem 20, we make the following 
conjecture: 

Conjecture 22. When A T ^ p\ there exists a catalyst state such that ® p\ P -5- T ® p M if and only if 

S t (A) >S t (n),Vte [l/2,oo). (52) 
It would also be desirable to determine the validity of Conjecture 19. 
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Appendix A: PROOF OF THEOREM 13 



The rank constraint 

rank X!''"^'' 1 X! ' '-'"'/ I ( A 1 

\i>j i>j 

can be rewritten as 

rank g ( Wj + (|ij)(ij| + |ji)0*l) / 2 + X A*<i|**X**l + X (^j - *«) (I**)0jI + l#X«l) / 2 = L ( A2 ) 

\i>3 i i>3 J 

The first sum and the remainder of the expression of the operator have disjoint support. The first sum can only have 
even rank, so to satisfy the condition we must have t i} - — —fiij (so that it's rank is zero). The remaining terms are 
then J2i ^u\ii)(ii\ + Z)»>j Mij (N)O'j'l + \ which has rank one if and only if fx l3 — UiUj, for some u e K d . 

Making the change of variables Xi = \/XiUi, the rank-constrained version of the optimisation problem (28) is 
therefore equivalent to 

JC2S1/2W _ 1 

Ti (K; X) - + max{A (a:) \x € S} (A3) 

K z — 1 

where 

A (-) = E|^-E|fT (A4) 

i>3 i>3 

and S is the hypercuboid defined by 

\xi\ < \/% for all i e {1, (A5) 




Figure 2: Illustration of the notations Fc and .He- 
rn the following we shall use some terminology from convex geometry ('face', 'supporting hyperplane'), for defi- 
nitions see Brondsted [14], for example. The function A is differentiable, so its maximum in the set S, A*, is attained 
by either an extremal point of S or a stationary point of the function on a face of S. 

If one of the components of x is negative, then changing its sign does not decrease A (x), so A* must be attained 
by a point on one of the faces of S in the set {F c \C C {1, d}}, where F c is the intersection of S with the supporting 
hyperplane H c defined by x t = +V^i, i e C (see Fig. 2). 
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The restriction of the function A to domain H c has a single stationary point at 



= / (K + \C\ - d) for all j e {1, d} \ C, 



(A6) 



where the value of A is equal to 



K 



{K 2 - 1) (K + c - d) 
The restriction of A to face F c contains a stationary point if and only if 



\\iec I iec J 



(A7) 



fev^j / (if + |C| - d) <min{ v / A~|*e{l,...,d}\C}. 



(A8) 



Lemma 23. Suppose we have two subsets of {1, d}, Ci and C 2 , both of size c, which differ by one element such that C\ = 
C U n\ and C 2 — C U n 2 . If A on Fc ± has a stationary point and X ni > A„ 2 , then A on F c , 2 also has a stationary point and 
A C2 > A Cl . 

Proof. Part (i) is trivial. For part (ii), since n 2 is not in C\, the fact that A on F Cl has a stationary point implies that 



J2 V^i < ( K + c - d ) V^r. 
iec 



and therefore 



(K + c-d) (K 2 -l)(A C2 -A Cl )/K 



= (K + c-d)(X ni -A n J- ^ V / V + - (^fiZ 

= (y/K~! - v 7 ^) I (K + c-d- I) (Vv + v^) -2 ^ 
> (a/V/V^- 1 ) (>/V + V^) >0. 



(A9) 

(A10) 
(All) 

(A12) 

(A13) 

□ 



If we assume, without loss of generality that A, = Aj, then as result of this lemma, if A* is attained by a point 
in the relative interior of an m-dimensional face then it must also be attained by the m-dimensional face with C — 
{1, d — m}. Therefore the A* is equal to Ac*, where C* = {1, c*} with c* being the smallest value (giving the 
largest dimensional face) such that A has a stationary point on F C ' or if no face has a stationary point then c* = d 
and the maximum occurs at the extremal point x = (y/X~i, The result follows. 



